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We consider the effect of adding quantum dynamics to a classical topological spin liquid, with 
particular view to how best to detect its presence in experiment. For the Coulomb phase of spin 
ice, we find quantum effects to be most visible in the gauge-charged monopole excitations. In the 
presence of weak dilution with nonmagnetic ions we find a particularly crisp phenomenon, namely 
the emergence of hydrogenic excited states in which a magnetic monopole is bound to a vacancy at 
various distances. Via a mapping to an analytically tractable single particle problem on the Bethe 
lattice, we obtain an approximate expression for the dynamic neutron scattering structure factor. 


The quest for spin liquids is an important enterprise in 
strongly correlated many body physics in an era when 
a huge amount of theoretical interest has focused on 
forms of order outside the canonical broken symmetry 
paradigm m- The search involves identifying relatively 
simple Hamiltonians that host spin liquids and finding 
experimental systems and signatures—the latter being 
more elusive than in Landau ordered systems. Indeed, at 
this point the list of experimental systems where there is 
strong evidence of spin liquid behavior is small. Among 
them is the celebrated spin ice system, arising in some 
rare earth pyrochlore magnets, which exhibits a U(l) 
spin liquid and excitations which are condensed matter 
analogs of Coulombically interacting magnetic monopoles 
[5]. Spin ice is truly special at this point in hosting a 
three dimensional spin liquid, but, owing to large mag¬ 
netic moments, is limited to the classical regime, in which 
coherent quantum dynamics appears to play little role. 

Logically, much recent interest has focused on looking 
for quantum generalizations of spin ice. There are several 
candidate materials for quantum spin ice behavior, such 

as Tb 2 Ti 2 07 [HI1171 . Yb 2 Ti 2 07 [10H12] and Pr 2 Zr 2 07 m, 

but an unambiguous experimental signature of quantum 
spin ice has been lacking. Logically, much recent theo¬ 
retical work has focused on looking for quantum gener¬ 
alizations of spin ice in which quantum fluctuations can 
lead to a fully quantum U(l) spin liquid [2H2D] • 

Here we investigate the addition of quantum fluctua¬ 
tions to spin ice but in a different limit which is, plau¬ 
sibly, of relevance to existing materials. Fundamentally, 
we wish to understand the leading order effects of adding 
quantum dynamics about the classical spin ice limit. As 
we will detail below, this has a parametrically larger ef¬ 
fect on monopole motion than on monopole-free ground 
states so the leading manifestations of quantum fluctua¬ 
tions appear when monopoles are present. 

We begin this program by studying the simplest man¬ 
ifestation of the quantum mechanics of monopoles—a 
striking effect that appears in the response of quantum 
spin ice to the introduction of a vacancy or missing spin. 
We find that the lowest lying excited states in the vicinity 
of the vacancy resemble those of hydrogen modulo lat¬ 


tice induced mixing—they involve a magnetic monopole 
bound to the impurity site into an infinite set of levels. 
In the presence of a dilute set of such impurities, these 
states give rise to a characteristic signature in neutron 
scattering at low temperatures which we discuss. Read¬ 
ers may note the family resemblance of these hydrogenic 
monopole states to hydrogenic states in doped semicon¬ 
ductors m and references therein], although we caution 
that the details have crucial differences. We also note 
that the response of spin liquids to impurities is of broad 
interest as a diagnostic of their internal dynamics: what 
happens when you dope a spin liquid is the fundamental 
- and to date largely unresolved - question of the RVB 
theory of high temperature superconductivity [25]. 

In the balance of the paper we begin by briefly re¬ 
viewing how the dynamics of quantum spin ice can be 
formulated as the quantum mechanics of monopoles. We 
then concentrate our attention on the problem of a va¬ 
cancy spin and describe how it can be mapped to a good 
approximation to a monopole on a Bethe lattice interact¬ 
ing with a fixed Coulombic charge. This model leads to a 
family of hydrogenic bound states of the monopole along 
with a continuum band. In the technical heart of the pa¬ 
per we solve this problem and obtain an exact closed form 
solution for the onsite Green’s functions. We use these 
results to obtain the signature of the hydrogenic states 
in the structure factor of spin ice containing a dilute set 
of vacancy spins. We conclude with some comments and 
pointers to future work. 

Quantum dipolar spin ice: Our model Hamiltonian 

HqbSI = Hdsi + E t-S, (1) 


consists, firstly, of the classical dipolar spin ice Hamil¬ 
tonian, defined for Ising spins Si living on the sites of 
pyrochlore lattice and pointing along the local easy axis 
joining centers of neighboring tetrahedra [23] [23] and ref¬ 
erences therein]: 
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FIG. 1. Spin ice projected onto a plane, with each vertex 
of the resulting square lattice in (a) corresponding to a py- 
rochlore tetrahedron, (b) A missing spin gives rise to a +1 
and a —1 charges; (c) flipping one of the majority spins ad¬ 
jacent to the vacancy creates a bulk charge +2, inverting the 
sign of one of the vacancy charges; (d) the bulk charge prop¬ 
agates in the system through further spin flips, while the net 
charge of vacancy is approximated as a single —2 monopole. 


ented perpendicular to the local easy axis, which adds 
the simplest quantum dynamics in the form of single spin 
flips. This simple form is convenient for a first theoretical 
analysis, for a more complete symmetry-based analysis of 
quantum terms in the Hamiltonian, see pan;]. 

Ghost spins and the Bethe lattice: At this point we 
switch from a spin description to that referred to as the 
dumbbell model [3] which we quickly review. Each spin 
is replaced by a pair of magnetic charges ±q m = /i /ad 
of opposite sign, where ad is a vector pointing between 
the centres of neighboring tetrahedra. By summing up 
the net charge at the center of each tetrahedron ( Q a = 
IlieaQi = 0, ±2q m , ±4g m ), we can replace the dipolar 
piece of the spin Hamiltonian © by 
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Coulomb interactions between charges on the diamond 
lattice, with Vo = 4 % \ f 1 + the cost of creating a 

monopole, which can be shifted by a nearest-neighbor 
exchange term. In spin ice ground states, Q a = 0. 
Flipping a spin in a ground states yields a pair of mag¬ 
netic monopoles of charges ±2g m on adjoining tetrahedra 
which can then move apart via further spin flips at finite 
cost in energy. Each charge has three majority spins that 
are all pointing in or out of the tetrahedron, and a single 
minority spin pointing in the opposite direction. 

With this in hand let us discuss the energetics of sub¬ 
stituting a magnetic ion on the pyrochlore lattice by a 
non-magnetic impurity, see Fig. [T] Removing a spin 


from a classical spin ice state (Fig. leaves behind 
two monopoles of charges ±q m [26] (Fig. So- A bulk 
monopole with charge ±2 q m can be ‘emitted’ by the va¬ 
cancy via flipping one of the two majority spins at each of 
the tetrahedra adjacent to the ghost spin (Fig. &)• The 
bulk monopole is then free to move around in the system. 

We define a quantity that can be thought of as the ioniza- 

2 2 

tion energy of the vacancy: I = ^^ 2 + 2uo^. This is the 

total energy cost of an emitted monopole moved out to 

infinity. Emitting another monopole into the bulk would 

2 

cost additional energy of the order of 2i>o ^ j so the lowest 

a d 

energy charged excitation in the presence of an isolated 
vacancy is a single monopole with charge ±2 q m . Once 
the monopole is emitted into the bulk, it is free to hop 
through fluctuation-induced spin flips, while a net charge 
of opposite sign ^2q m remains at the vacancy (Fig. ^Ji). 

Adding quantum dynamics via the transverse field in 
Eq. ([I]) has only a weak effect on the ground states, as 
connecting two of them requires flipping spins in closed 
loops, minimally six of them on a hexagon of the py¬ 
rochlore lattice. Near the classical limit, vq t, such 
processes come with a prohibitively small energy scale, 
~ t 6 /v g. By contrast, for a state containing a monopole, 
the lowest order effect - a monopole hopping onto a 
neighboring tetrahedron by flipping a majority spin is 
parametrically stronger: linear in t\ 

Thus, in experiment, the most promising place to 
see quantum effects in spin ice is in the gauge-charged 
monopole excitations, rather than its gauge-neutral gap¬ 
less emergent photons. Analogous considerations apply 
in the proximity of a vacancy, where we focus on the case 
of a monopole emitted into the bulk (Fig. [lji), also with 
low-order signatures. For this reason, here we perform a 
quantum calculation for the monopole states, and do a 
thermal sum over the nearly degenerate spin ice configu¬ 
rations. 

We treat the problem as that of two Coulombic 
charges, one of which is stationary. As the charge prop¬ 
agates through the bulk, it changes the spin ice back¬ 
ground. This process is difficult to capture exactly, but 
fortunately it is possible to make considerable progress 
via an effective model that we describe next. From this, 
we are able to extract the bound states in considerable 
detail, followed by a continuum band, much as we would 
expect for the Hydrogen atom. 

In order to investigate the problem of an isolated va¬ 
cancy that has emitted a free monopole into the bulk, 
we switch to the state lattice description m- First, con¬ 
sider a new basis of the following (classical) states: a 
spin ice state with a vacancy, which we label 10), and 
states with an emitted monopole in the bulk, connected 
to |0) through single spin flips. Next, each site of the 
state lattice represents one of the basis states |n); while 
bonds connect those sites whose corresponding states are 
connected by single spin flips. Apart from site 0 (repre- 
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FIG. 2. The Bethe lattice describing state space. Its root, 
n = 0, corresponds to the unionized vacancy state |0). 


senting |0)), the state lattice is trivalent. It can be shown 
that the smallest closed cycle in the state lattice of disor¬ 
dered pyrochlore spin ice has length 20 [25]. We therefore 
approximate the state lattice by a cycle-free infinite Cay¬ 
ley tree (the Bethe lattice) rooted at site 0 (Fig. [2]). The 
monopole propagating in real space corresponds to a sin¬ 
gle particle hopping on this lattice in the presence of the 
Coulomb potential: 

#i°> = -*'52 l TO ) ; H ^ n '> = ( I+ J~) l n M 52 l 77 ^ 

m= 1 ' ' m—1 

(4) 

where the sums run over states reached by flipping ma¬ 
jority spins of the monopole. d n denotes the distance 
between the vacancy and the monopole in the bulk in 
units of ad, such that C/d n is the attractive Coulomb 
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potential [C = between the two charges. In 

conventional spin ice, the cost of having a monopole is 
larger than the magnitude of the Coulomb interaction 
between two charges, I > \C\. For concreteness, we use 
C = — J/3; t = //10 in the following. Since we restrict 
ourselves to a particular starting spin ice configuration 
and omit other degenerate ice states from the discussion, 
the mapping from Eq. |l]) to Q is accurate to 0(t 5 /v g). 
Our final approximation concerns the distance between 
monopoles. Since the four sites at the first generation of 
the Bethe lattice correspond to the bulk monopole being 
one spin flip away from the vacancy, it is natural to ap¬ 
proximate d n in Eq. Q by the generation of the Bethe 
lattice n. This definition fails to be exact already beyond 
0{t 2 ), but should work sufficiently well in the /, |Cj t 
regime, when the bulk monopole prefers not to move too 
far. In return for these approximations, we are able to 
solve exactly our idealized model, that of a single par¬ 
ticle hopping on the Bethe lattice in the presence of a 
Coulomb potential / + C/n for n > 0. 

The Bethe lattice problem: We calculate the diago¬ 
nal elements Gu(tx>) of the lattice Green’s function to infi¬ 
nite order in t [291130] . We find [28] that each Gu{u) can 
be written down in terms of a finite number of G^(co), 
infinite sums involving particle hopping from a site at 
generation k to sites at generations g > k. The latter 
have a closed form expression in terms of the Gauss hy¬ 
pergeometric functions F 2 (a, b , c, z) [31]: 
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where x 2 = — This yields the exact expression for 
any of the diagonal elements of the Green’s function; for 
instance at the root site 

SooM = (w — (w — -0) -1 • 

The full Green’s function yields the energy levels via its 
poles and the local densities of states for each Bethe lat¬ 
tice generation, proportional to its imaginary part. The 
local density of states at site 0 in Fig. [3] indicates that 
indeed there are bound states followed by the contin¬ 
uum energy band. While the classical ground state (a 
spin ice state with a vacancy) would have zero energy, 
the ground state energy of the quantum problem ojq is 
lowered due to the hopping t. Low-lying excited states 
are separated from the ground state by a gap, which is 
also decreased from the classical value I through hop¬ 
ping and Coulomb attraction. They accumulate below 
the edge of the continuum band, located at L = I—V&t 2 - 
In the Bethe lattice problem, the band of the extended 


states, of width linear in t, is confined to the region 
I — \/8 k 2 < u> < I + V&t 2 . (Introducing closed cycles 
into the lattice has the effect of adding band tails, ex¬ 
tending beyond these edges.) 

Signatures of monopoles in neutron scattering: 

One of our central results is the dynamic structure factor, 
defined as 

S(Q, Au,) =525(E f -E i -Au,)\J2(f\Sp)e i *«\ 2 . (6) 

/ R 

In order to extract the information that is most relevant 
to spin ice experiments from the Bethe lattice model, 
we calculate a one dimensional version of Eq. (|6|, aver¬ 
aged over all directions of q. Such a quantity, S(q,Aui), 
can be measured directly in a powder averaged neutron 
scattering experiment. The details of our calculation, 
carried out in the limit of dilute nonmagnetic impuri¬ 
ties, are given in the Supplemental Material [25], The 
dynamic structure factor S(q, Aw), plotted in Fig. [4[a), 
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FIG. 3. Local density of states at site 0 of the Coulomb prob¬ 
lem on the Bethe lattice, for C = —1/3', t = 7/10. Bound 
states (red) appear as sharp peaks; the lower edge of the con¬ 
tinuum (blue) is labeled L, the classical ground state energy 
by 0 and the ionization energy by 7. 



FIG. 4. Left: dynamic structure factor S(q,Auj) for pow¬ 
der averaged neutron scattering. Each line is multiplied 
by Exp[0.7|n — 1|], where n labels the all-even energy lev¬ 
els with n = 0 the ground state. Right: Lineshapes 
S(q, Aui n )/S(0, Auj n ) for n = 0,1, 5,10. C = —7/3; t = 7/10 
throughout. 


has sharp features signaling the presence of bound states. 
The structure of the lines gives direct information about 
the character of the ground and excited states. The most 
visible signatures show up in elastic scattering and at the 
energy transfer equal to the difference between the first 
excited state and the ground state. For a well-localized 
ground state, the matrix elements between S'ili) and ex- 
cited states at higher energies (bound to the vacancy at 
distant radii) give rise to peaks whose structure is essen¬ 
tially identical up to a scale factor, as shown in Fig. S b )- 
Note that for 7 = 0, the signals corresponding to n ^ 1 
would be absent, vanishing as powers of t. Their presence 
thus yields direct evidence of the existence of quantum 
dynamics. 

Conclusion and outlook: We have studied in detail 


the properties of magnetic monopoles in dipolar quantum 
spin ice. We have demonstrated that these are the prime 
indicators of the presence of quantum dynamics. In the 
presence of nonmagnetic impurities we have found both 
sharp liydrogenic bound states as well as a broad contin¬ 
uum energy band. While we believe these results to be 
robust, there is clearly much scope for further, presum¬ 
ably numerical, modeling taking into account the detailed 
lattice structure, as well as any material specific single¬ 
ion physics and terms in the quantum Hamiltonian. 

The quantum dynamics of a pair of monopoles presents 
a more difficult problem due to the pair’s center of mass 
motion. We are planning to address this issue, as well 
as clarify the detailed character of the continuum band 
of states in the vacancy problem, in future work. Addi¬ 
tionally, despite neutron scattering being the method of 
choice for investigating magnetic materials, local disorder 
is an attractive subject for other types of experimental 
probes, such as nuclear magnetic resonance. While such 
techniques are beyond the scope of this work, our the¬ 
oretical model can also be employed for calculating real 
space quantities accessible by the local probes. 
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SUPPLEMENTARY MATERIAL 
Minimal closed cycle in the state graph of pyrochlore spin ice 



FIG. 5. A fragment of the real-space diamond lattice that the magnetic monopoles of pyrochlore spin ice reside on. Spins are 
depicted as arrows pointing along the lattice edges. Left: the shortest closed cycle in the state graph involves flipping spins 
along two hexagons that share two edges in common (red). Right: an example of a spin configuration, that gives rise to the 
shortest possible cycle on the state lattice, with two flippable hexagons and a magnetic monopole at one of the sites where they 
intersect. 

One of the arguments for using the Bethe lattice to approximate the state graph of pyrochlore spin ice is that the 
closed cycles present in the actual state lattice are long and thus can be neglected. The spins can be thought of as 
pointing along the edges of the diamond lattice [5j. Its smallest closed loop has length six, so inverting spins along a 
flippable hexagon (such that six spins are arranged to point head-to-tail) brings the charge back to where it started. 
The spin configuration is changed by such a process, and one needs to flip the same six spins a second time to bring 
the system back to the initial state. However, this process constitutes a self-retracing path on the state graph, rather 
than a closed cycle. 

Introducing a loop into the state graph requires adding another hexagon into the real-space picture. The shortest 
possible loop is obtained by maximizing the overlap between the two hexagons: they share two edges as shown 
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in Fig. [5] A non-trivial closed loop on the state graph involves ten spins on two such overlapping hexagons, and 
positioning a magnetic monopole on a site from which it can hop along two different paths [charge —2 q m on site 2 
in Fig. [ 5 ]. Each spin can be labeled by two numbers that correspond to the diamond lattice sites at the two ends of 
the dumbbell. An example of a closed cycle in the state graph, written down as a sequence of spin flips as labeled in 
Fig- [5} is then: 


(23) (34) (45) (56) (61) (12) (29) (98) (87) (76) (65) (54) (43) (32) (21) (16) (67) (78) (89) (92). 

This closed cycle has length 20. We would like to emphasize that this is the shortest possible loop in the state lattice 
of disordered classical ice states. Other spin configurations have closed cycles that either have the same length, or are 
even longer. 


Lattice Green’s function: the continued fraction method 


In order to solve the Coulomb problem on the Bethe lattice, we employ the continued fraction approach [29] [30]. 
The strategy is to start with a perturbative treatment (where H 0 is the potential, and nearest neighbor hopping t is 
the perturbation) and proceed to solve the problem exactly by calculating the self-energy to infinite order in t in the 
particle’s lattice Green’s function. Recursively solving Dyson’s equation 

Q(uj) = Gq{oj) + Go(ui)Y,(u))G(uj) 

(where E(w) is the self energy) gives rise to continued fractions for the Green’s function. For instance, the diagonal 
element of the Green’s function for site 0 is: 


£oM = 


w - S 0 (w) 


4 1 2 


u> — 


■-I-C- 


2 1 2 




The self energy Eo(w) above is given by the sum of all paths on the lattice going away from site 0 and back to it. 
The elements of the Green’s function for sites at all generations can be defined in terms of a finite number of Gl f (w), 
infinite continued fractions involving hops only to generations higher than k: 


SfM = 
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Using a theorem by Ramanujan }3l[ . we arrive at the following exact closed form expression for the continued fraction 
above: 


SiTM = 


2 k/u 
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where x 2 = — and F 2 (a,b,c,z) is the Gauss hypergeometric function. 


Calculation of the dynamic structure factor 

Quantum numbers of the Bethe lattice eigenstates 

Consider a set of only three sites: site 1 at generation n and sites 2 and 3, connected to it, at generation (n + 1). 
Exchanging sites 2 and 3 leaves the Hamiltonian @ invariant. Now examine what happens when Eq. 0 acts on 
two states: |’F,s) and |SR^i), symmetric and antisymmetric combinations of the particle being at sites 2 and 3. The 


























7 


P(n) 



FIG. 6. Probabilities of the four lowest energy states in the all even sector, plotted as a function of the Bethe lattice generation. 
C= -7/3; t = 7/10. 


diagonal part of the Hamiltonian will be the same for both states, and both states allow for the particle to hop to 
higher generation sites. However, hopping back to site 1 is eliminated by I'I'a)) since the contributions from sites 
2 and 3 cancel out. Similarly, each Bethe lattice generation n > 1 can be assigned an even/odd quantum number, 
corresponding to the symmetries that involve exchanging left and right sites at each generation. If a particle starts 
at a state that is odd at k th generation, it can only hop to sites at generations n > k. It follows that the particle’s 
wavefunction has zero amplitude at all generations n < k. Therefore, only the all even states have a nonzero amplitude 
at the origin. 


Powder averaged dynamic structure factor 

The dynamic structure factor is defined as 

S(q, Aco) = $>(£/ ~ E i~ A w )|^r</|S+|*}e^| 2 . 

f R 

We calculate the exact expression for S(q. , A lj) that is averaged over all directions of q for the Bethe lattice problem. 
This gives us the approximate dynamic structure factor S(q, A lj) that can be measured in a powder averaged neutron 
scattering experiment on quantum spin ice, in the limit of weak dilution. 

The ground state |i) of the Coulomb problem on the Bethe lattice can be written down in the following form: 

I*) = «|0) +/3|1) +7l 2 ) + • ■ • 

where the Greek letters represent the ground state’s amplitude per lattice generation. In the course of a scattering 
experiment, a neutron coming in flips a spin. Since we restrict ourselves to having at most one free monopole with 
charge ±2 q m , the spin flip either nucleates a monopole at the vacancy, eliminates it, or shifts an existing monopole 
by a lattice spacing: 

J2 E<rii S±\i) = a (4e i<? 11)) + /3 (e* 9 |0) + 2e i2g \2)) + 7 (e i2q \\) + 2e* 3 ? |3)) + ... 

R 

Since the ground state belongs to the all even sector, the signs of the amplitudes for the sites at the same generation 
are equal, and the contributions from the spin flips add up (hence the numerical prefactors in the expression above). 
When we take the inner product of YIr e l ^' R Si\i) with excited states which contain odd quantum numbers, however, 
the hops to sites at the same generation cancel out. Therefore, only the all even states contribute to the powder 
averaged structure factor. We can see how these states are distributed in real space by plotting their probabilities as 
a function of Bethe lattice generations n (Fig. [6]). In particular, the ground state is localized near the origin, which 
leads to the scattering intensity being high for lower energy states, and low for the states that accumulate near the 
continuum band edge, whose wavefunctions near the vacancy differ in amplitude, but not much in shape. 






